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Abstract 

We present a lattice study of the phase transitions at zero and nonzero temperature for the SU{3) gauge theory with 
a varying number of flavours in the fundamental representation of the gauge group. We show that all results are 
consistent with a lower edge of the conformal window between Nf = 8 and Nf =6. A lower edge in this interval is in 
remarkable agreement with perturbation theory and recent large-A arguments. 
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1. Introduction 

Physicists are familiar with two paradigms for confor¬ 
mal symmetry loss in four-dimensional interacting quan¬ 
tum field theories: Quantum Electrodynamics (QED), whi¬ 
ch is infrared free, and Quantum Chromodynamics (QCD), 
which loses conformal symmetry in a highly non trivial 
way, leading to two manifestations of one single breaking 
phenomenon: asymptotic freedom and confinement 

A richer dynamics can be realised in the presence of 
non-trivial, i.e., interacting, ultraviolet or infrared fixed 
points of the quantum theory, and may turn out to play a 
role in unifying the standard model of particle physics with 
gravity. The possibility that such a scenario is realised 
somewhere between the electroweak symmetry breaking 
scale and the Planck scale has recently attained strong 
theoretical and experimental appeal. The search for re¬ 
normalization-group fixed points in candidate theories for 
particle physics and cosmology also shares a wide range 
of tools and aspects with the study of quantum critical 
phenomena in condensed matter systems, in three as well 
as lower spatial dimensions. 

The simplest step beyond the QCD paradigm is the 
so-called conformal window: a family of zero-temperature 
theories in the phase diagram of non-Abelian gauge theo¬ 
ries ranging over a number of massless flavours N'j <Nf < 
Nuv 

in a given representation of the gauge group. The 
critical number of flavours Nj marks its lower edge, and 
the upper edge where ultraviolet freedom is lost. The 
beta-functions f3(g) of theories inside the conformal win¬ 
dow have a zero at some coupling g = g*, where the theory 
has a non-trivial infrared fixed point (IRFP) [TJ [5] and it 
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is conformal; they are negative for 0 < g < g*, and zero 
at 5 = 0 (ultraviolet freedom). Above the conformal win¬ 
dow, f3(g) > 0 for g > 0 implies infrared freedom opening 
up the interesting possibility of asymptotic safety, i.e., a 
non-trivial ultraviolet fixed point (UVFP) [3]. 

Many questions related to non-trivial fixed points and 
the emergence of the conformal window — its location in 
the parameter space and accompanying signatures — are 
genuinely nonperturbative, and require nonperturbative 
strategies. Gauge/gravity duality [1] can come to rescue 
whenever the gauge theory is conformal, as in the confor¬ 
mal window, or near-conformal and still deconfined, as in 
the high-temperature quark-gluon-plasma. This seems to 
be not true in the confining and asymptotically free phase 
of QCD [5 that precedes the conformal window; the re¬ 
cently proposed solution for the scalar glueball corre¬ 
lator in the ’t Hooft limit of large-iV QCD tells that its 
momentum dependent logarithms cannot be reproduced 
by AdS/CFT realisations so far. 

In this letter we constrain the lower edge of the confor¬ 
mal window for the SU(3) gauge theory with Nf flavours 
in the fundamental representation (i.e., many-flavour mass¬ 
less QCD) by means of a lattice study, a genuinely non¬ 
perturbative approach, tailored to the physics problem at 
hand. Knowing the location of the lower edge is one way 
to establish how far large-fV predictions, or perturbation 
theory to a given loop order are from the complete theory. 
It is also essential in order to correctly identify the proper¬ 
ties of many phenomenological models for particle physics 
with a composite spectrum. 

2. The lower edge of the conformal window 

Figure summarises the results of this work and ex¬ 
plains its strategy: theories with Nf < Nf have sponta¬ 
neously broken chiral symmetry below a critical temper¬ 
ature Tc > 0, above which it is restored. Inside the con¬ 
formal window, theories have exact chiral symmetry and 
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Figure 1: Lattice phase diagram Nf — 1//3, /3 = 10lg^, the lattice 
coupling. Data for Nf = 8, ..., 12 (red circles) are midpoints of sharp 
zero-temperature crossovers in Fig. Best fit to data (solid red) 
identifies the linear in 0 bulk line, extrapolated to Ng < 8 (dotted 
grey). Nf = 7 crossover in Fig. plfalls on the line (grey circle). Bound 
on N'j from this work is sketcned (grey band). Nj = 4, 6 thermal 
transitions occur away from the bulk line, with volumes 24^ X Nt, 
Nt = 6 (orange), 8 (green), and 12 (blue), right to left. 

are deconfined at all temperatures. Hence, one way to 
determine Nj is to follow the chiral symmetry breaking 
pattern, with varying Nf and temperature, for the theory 
regularised on a Euclidean spacetime lattice. For Nf< Nf, 
a thermal chiral symmetry restoring phase transition oc¬ 
curs at some Tj, = l/(a(g£)A^t) > 0, with a{g1) the lattice 
spacing at the critical lattice bare coupling gl on a lat¬ 
tice volume Nj^ x Nt, with temporal extent Nt<^Ni. For 
Nf ^ Np the lattice theory exhibits exact chiral symme¬ 
try at all temperatures, i.e., any Nt, and for all values of 
gL in the interval 0 ^ gL ^ 9 l^ with ^ sufficiently large 
coupling where chiral symmetry will eventually be broken. 
This is the bulk phase transition. 

The red solid line in Fig. [^is the line of such transitions 
for varying Nf. In any given renormalisation scheme, it 
manifests the fact that fermion screening is increasingly 
effective for increasing Nf inside the conformal window; 
this Nf dependence is a leading order effect separating 
two phases with different underlying symmetries, different 
in nature from any lattice artefact that could occur for 
Nf < Nf inside a single chirally broken phase on coarse 
lattices. At finite lattice spacing a, the bulk line can be 
seen as the Nt = l/(oT) —>■ oo limit of an A"t-finite family 
of chiral phase transitions that exists for all Nf > 0. The 
line therefore stops at iVji, because no T = 0 chiral phase 
transition occurs for Nf <Nf. 

With this paradigm, we have studied the Nt depen¬ 
dence of chiral phase transitions for 4^ Nf ^ 12 in QCD. 
The setup of the simulations is the one of [5]. Through¬ 
out, we have not observed an anomalous behaviour that 
could hint at consequences of the fourth root of the fermion 
determinant for staggered lattice fermions. We have stud¬ 
ied the chiral condensate (ipilj), order parameter of spon¬ 
taneous chiral symmetry breaking, the disconnected chi¬ 
ral susceptibility Xdisc, and the connected chiral cumu- 

lant i?^=Xconn /X7r = (9(l^^/>)/l9mvalence)/(('*/''0)M)- T^is 

quantity is a powerful probe of chiral symmetry: in the 
chirally broken phase, the pseudoscalar lowest-lying state 


Figure 2: Sharp crossovers of the chiral condensate (lattice units) 
for Nf = 12 [ 7 ] to 7 (left to right). Nf = 11,10, 9 at 12^ x 24, Nf = 8 
at 16^x32 (circles, solid line), 24^x12 (triangles, dotted line), 24^x6 
(squares, dashed line), Nf = 7 at 16^x32; all with bare lattice mass 
am = 0.01. Inset shows data for two branches of possible hysteresis 
loop for Nf = 7. Edge of exotic phase displayed for Nf = 12, 9, and 
Nf=8. 

is a Goldstone boson and its vanishing mass in the chi¬ 
ral limit, together with the non-degeneracy of chiral part¬ 
ners, guarantees that —>■ 0 in the chiral limit; in the 

chirally restored phase, the degeneracy of the scalar and 
pseudoscalar chiral partners implies that > 1 in the 
chiral limit. We have also studied the Polyakov loop L 
at finite temperature. Despite it not being a true order 
parameter in the presence of fundamental fermions, it is 
observed to typically retain the correct features related to 
a deconfined or a confined phase. We used volumes with 
Nt < Ni and aspect ratio NfNt ^ 2, as well as zero tem¬ 
perature volumes with Nt = 2Ni , at a fixed bare fermion 
mass 0771 = 0.01. 

Remarkably, the sharp crossovers in Fig. [^for Nf = 12 
to Nf = 8 are Nt independent. We do observe a reduced 
crossover for Nf = 7, an almost closing gap, and signs 
of hysteresis in the interval /3 = [4.275,4.3] right on the 
bulk line of Fig. We defer the fate of Nf = 7 for a 
future study. We have also extrapolated the bulk line to 
Nf = 16 obtaining /3c ~ 1.2, consistently with previous 
studies [7] carried out with different lattice actions and 
heavier fermions. 

A genuinely different behaviour is observed for the Nf = 
4, 6 theories, whose markedly W-dependent chiral symme¬ 
try breaking crossovers for varying Nt < Ni are reported 
in Fig. Already at Nt = 6, the crossover occurs at a 
coupling weaker than the one predicted by the linear in 
/3 extrapolation of the bulk line. It moves to weaker cou¬ 
pling for increasing Nt, consistently with a thermal nature 
of the transition. At the same time, no zero-temperature 
{Nt ^ Ni) phase transition or analytic crossover is observed 
to occur along the extrapolated bulk line, nor at weaker 
coupling. 

These results provide evidence that the lower edge of 
the conformal window lies between Nf=8 and Nf = 6, the 
main result of this work. In what follows, we further scru¬ 
tinise the Nf = 8 and Nf=6 theories, discuss implications 
and future directions. 
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Figure 3: Nf=S observables (lattice units), top to bottom: {'ip'ilj), 
connected i? 7 r, Xdisc? (ReL), 24^ x = 6 (red squares), 

A^t = 12 (green circles), and 16^x32 (blue triangles), mass am = 0.01. 
{ReL), Nt = 12 (HYP-smeared) rescaled for better visibility. Vertical 
lines indicate the location of the bulk transition as in Fig. ^ and 
the edge of the exotic phase for Nt = Q, 12, 32, left to right [We have 
obtained identical results with volumes 16^ x32 and 12^ x24, which 
suggests that finite volume effects are under control for Nf = 8, as 
well as Nf > 8]. 


3, The Nf — 8 theory and the exotic phase 

Figure shows strikingly the Nt independence of the 
strong coupling crossovers of o-nd their perfect over¬ 
lap with those of R-jr-, and with the A^^-independent peaks 
of Xdisc- These signals are consistent with exact chiral 
symmetry for the Nf = 8 theory at all temperatures, and 
the occurrence of a chiral symmetry breaking bulk tran¬ 
sition at sufficiently strong coupling, as reported in Fig. 

For A^y = 12, an extensive analysis for varying volumes 
and fermion masses has suggested a first order nature of 
the transition also away from the chiral limit, for a bare 
fermion mass am = 0.01 nm. 

Towards weaker coupling, Fig. and Fig. show a 
sequence of Vj-dependent reduced crossovers of the chiral 
order parameter, coinciding with a dip of the connected 
Rt^ for all Nt^Ni, and accompanied by a sharp crossover 
in the finite temperature {Nt < Ni) Polyakov loop, which 
remains zero at stronger coupling. We note that this sec¬ 
ond crossover in (V't/') is clearly visible for Nt = 6, whilst 
its signal is suppressed at larger Nt = 12,32, due to fermion 
mass effects at weaker coupling for fixed am. Two facts 
are important. Firstly, signals survive the zero tempera¬ 
ture limit, and their Nt dependence ceases for Nt ^ 12; 
this excludes their thermal nature. Secondly, the absence 
of corresponding Xdisc peaks is consistent with restored 
chiral symmetry at all temperatures, and coupling smaller 
than the bulk coupling. These signatures are recognised 
to identify the weak coupling edge of the so-called exotic 


phase, a genuine lattice artefact extensively studied mm 
in the Nf = 12 theory, where it is observed to precede the 
bulk transition when improvement of the lattice fermion 
action is used [7]. They also suggest that the transition 
interpreted as a thermal transition for Nf = 8 in [iiiin]- 
with a different improvement of the fermion action — is 
instead the edge of the exotic phase, and that both Nf=8 
and Nf = 12 are inside the conformal window. 

That the Nf = 12 theory is inside the conformal window 
and has exact chiral symmetry at weak coupling has been 
corroborated by the vanishing of the inhnite volume chiral 
condensate in the chiral limit m and a thorough universal 
scaling study of the spectrum of the theory [5], including 
all lattice results available at the time [TH [T5] . The study 
of Nf = 12 in [7], see also [16], has also provided consis¬ 
tent evidence that the exotic phase is a chirally symmetric 
phase, where the U(l) axial symmetry is also effectively 
restored. 

One might wonder whether the improvement of the 
lattice fermion action could have spoiled the thermal be¬ 
haviour of a theory below the conformal window, and the 
exotic phase as we observe it could have arisen inside a 
chirally broken phase. The answer is no, according to 
the study of = 12 with and without improvement |7j; 
both cases indicate that chiral symmetry is not sponta¬ 
neously broken and no thermal phase transition occurs. 
Furthermore, the lattice spacing a* at which the exotic 
phase would emerge, and the lattice spacing Oc at which 
a thermal phase transition would occur, carry an a priori 
different Nf and Nt dependence. It would be highly acci¬ 
dental to have a* < Oc for any studied Nt and all Nf ^ 8, 
in order to remove traces of a thermal transition, whilst 
preserving the shape of fVt-fixed lines for all Nf^6. 

The exotic phase observed here has its seeds in the 
conformal window, which genuinely differs from QCD. It 
appears inside the chirally symmetric Coulomb phase of 
the lattice regularised theory (see Figure 2 of |8]), where 
no continuum limit can be taken and improvement — used 
to reduce lattice artefacts at weak coupling— is not even 
approximately justified, because the theory is no longer 
ultraviolet free. It has a truly exotic nature: it is chi¬ 
rally symmetric but behaves as confining, oppositely to, 
e.g., finite temperature QCD with adjoint fermions, where 
a deconfined but chirally exact phase is allowed and in¬ 
deed observed. We add that this and similarly originated 
exotic phases may acquire physical relevance in discrete 
ANNNI models and systems such as graphene in varying 
dimensions. A complete microscopic understanding of the 
exotic phase is, however, still lacking. 

4. The Nf = 6 theory and asymptotic scaling 

A study analogous to Nf = 8 reveals a different sym¬ 
metry pattern for the Nf=6 and Nf=4: theories. Figure 
shows the A'^-dependent signal of the chiral observables 
and the Polyakov loop for varying Nt < Ni for Nf = 6, 
and no exotic phase. All crossovers overlap with the Xdisc 
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Figure 4: Nf=Q observables (lattice units), top to bottom: 
connected i? 7 r, Xdisc? (ReL), 24^ x Nt, Nt = 6 (red squares), 

Nt = S (green circles), Nt = 12 (blue triangles) and 16^ x 32 (parse 
diamonds), mass am = 0.01. (ReL), Nt = S and Nt = 12 (HYP- 
smeared) rescaled for better visibility. Vertical lines indicate the 
locations of the thermal crossovers as in Fig. ^ Vt =6, 8, 12, left to 
right. Inset for Xdisc- enlarged Nt = 12 peak versus zero temperature 
fiat distribution. 

peaks. For a comparison, the = 6 signals for Nf =6 in 
Fig. 0 and Nf = 8 in Fig. are equally sharp, but with 
genuinely different characteristics. 

The thermal nature of the Nf=6 and Nf=4: crossovers 
is corroborated by the two-loop asymptotic scaling study 
in Fig. similar in spirit to |12j . In particular, we com¬ 
pare the lattice data to the two-loop scaling curve 



Figure 5: Lattice data with Nt = Q, 8, 12, Nf = 6 (left) and ^^ = 4 
(right) as a function of the £^-scheme improved coupling 
compared to the 2-loop scaling curve ( 5 '£;- 2 /) in equation 0 , const = 
12 * R{g^(Nt = 12)), and the lattice-distorted curve (p£:-2/LD) in 
equation const = 12 + /? 2 >D( 5 '^('^t = 12)), = 0.006 (Nf=6) and 

h = 0.000D[Nf =A). Error bands are scaling curves for the extremes 
of the error on g^(Nt = 12). 

the fact that the latter theory has a lower critical tem¬ 
perature, hence a larger critical lattice spacing and larger 
lattice artefacts, for the same Nt. However, if contribu¬ 
tions to asymptotic scaling beyond two loops are signifi¬ 
cant (more likely for Nf=6 than TV/=4) the h parameter 
in the two-loop scaling will be affected by higher loop ef¬ 
fects. A complete four-loop analysis in the MS scheme is 
desirable, but we currently lack the full conversion of the 
lattice bare coupling to the MS coupling for the given lat¬ 
tice action. In using the asymptotic scaling formula valid 
in the massless limit, we have assumed that the relative 
shift of the (pseudo) critical couplings due to the nonzero 
fermion mass with varying Nt is within present uncertain¬ 
ties. A more refined study is certainly worthwhile, follow¬ 
ing the strategy used for the precise determination of the 
QCD pseudocritical temperature ^0] . 


^^ 3 E)Nt = = const, (1) 

with R{gE){=aAE) the two-loop asymptotic scaling func¬ 
tion and the E-scheme improved coupling m 

fe" = (l/3)/(l-(P)Mc/3), (2) 


with {P)mc the Monte Carlo determined zero tempera¬ 
ture plaquette at the critical coupling. The constant in 
equation Q is fixed by the data at ^( = 12. 

The leading-order lattice-distorted two-loop scaling {gE- 
2TLD) [18], which corrects for lattice artefacts due to the 
finiteness of a, 


Rld(.9e) = R ^{9e) 


1-h 


R\9e) 


R^{g%{Nt = 12)) \ 


(3) 


is in good agreement with the data for h A larger 

h is needed when the tadpole improved m (lattice bare) 
coupling is used in equation (§, /i = 0.024 (Ti = 0.0275) and 
/i = 0.018 (/i = 0.022) for Nf = 6 and Nf=4, respectively. 
The increase of h from Nf = 4 to Nf = 6 is consistent with 


5. Discussion 

This study consistently suggests that the conformal 
window of QCD opens between Nf = 6 and Nf = 8 fun¬ 
damental flavours. The devised strategy makes use of 
the underlying symmetries of the SU{3) theory and their 
breaking pattern with varying Nf and temperature, and it 
should be preliminary to the lattice investigation of other 
properties. 

Consider the hadron spectrum; it is numerically chal¬ 
lenging to distinguish between the spectrum just above 
and just below the lower edge, in the presence of a nonzero 
fermion mass and at finite volume. For example, the fea¬ 
tures of the Nf = 8 spectrum studied in [21] can be ex¬ 
plained by an underlying exact chiral symmetry, in agree¬ 
ment with this work: the lowest-lying scalar meson is light 
because it needs to be degenerate with its chiral partner, 
the pseudoscalar, in the chiral limit. One should also ob¬ 
serve that inside the conformal window residual finite vol¬ 
ume effects behave differently from QCD, given that only 
a long-range spin-dependent Coulomb potential is present. 
In fact, the study of Nf = 12 [5] has shown that spin-1 
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channels can be more affected than spin-0 channels by the 
finite size of the box, leading to an enhancement of the 
TTip/m^ finite volume ratio not dissimilar to Nf =8 in |21j . 

A lower edge of the QCD conformal window around 
Nf = 7 is not far from the prediction of two-loop per¬ 
turbation theory, it is in remarkable agreement with the 
large-result Nf/N = 5/2, for N = 3 [55] and with four- 
loop perturbation theory [53] . It is also consistent with the 
perturbatively small value of the fermion mass anomalous 
dimension of the Nf = 12 theory mm. because Nf = 12 
is not close to the lower edge of the conformal window. 
A more accurate study of iV^ = 7 at a smaller fermion 
mass should allow to determine if this theory is conformal 
or confining, while any form of preconformal behaviour, 
if it exists, is likely to be manifest in the Nf = 6 theory. 
The latter is not far from real world QCD and it remains 
an instructive playground for standard model extensions, 
including dark matter. 
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